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1. Introduction 

The low energy dynamics of QCD can be quantitatively understood computing the universal low 
energy constants (LECs) of the chiral effective theory [JTJ] . QCD results obtained with simulations 
on a space-time lattice can be matched with the effective theory, if the range of quark masses and 
energy scales in the QCD computations, are such that the higher order corrections in the effective 
theory calculations are well under control. A particular attractive framework in this respect is the 
so-called £ -regime In this regime the power counting for masses, momenta and linear size of 
the volume, is such that higher order corrections are suppressed in comparison with the standard 
/^-regime. Thus this regime provides an alternative, but also complementary, way to determine 
the LECs in comparison with standard p-regime matchings. The £-regime is properly matched if 
the volume in the QCD simulations, and correspondingly in the chiral perturbation theory (#PT) 
computations, is larger than the confinement radius and if the Goldstone boson correlation length 
of the system is smaller than the inverse linear size of the volume. Typically the usage of Ginsparg- 
Wilson ^ (GW) fermions is preferred because GW fermions have a natural definition for the 
topological charge, thus they allow to study the £ -regime in a fixed topological sector |Q] for which 
also #PT formulae have been worked out [§J. Wilson-type fermions can still, in principle, probe 
the £ -regime, sampling all topological sectors ^]. An important point is the interplay between 
low quark masses and a finite lattice spacing. It can be addressed using the concept of generalized 
chiral expansions applied to the £-regime [||, A second issue is the algorithm which ought to 
be used to practically simulate in the £-regime with Wilson-type fermions. In these proceedings 
we cover these two topics, and present preliminary results for the chiral condensate 1 £ and the 
pseudoscalar decay constant F. Additionally we discuss the systematic uncertainties which can 
affect the determination of these LECs in the £-regime. 

2. e expansion with Wilson fermions 

The £-regime has been introduced to cure the appearance of infrared divergences in #PT when 
the mass squared of the Goldstone boson (M n ) becomes smaller than the inverse size of the box 
(V = L? x T). The infrared divergences appear because the zero modes of the Goldstone bosons are 
treated on the same footing as the non-zero modes. To cure this problem Gasser and Leutwyler [Q] 
proposed to change the power counting (from /j-regime to £-regime) in this particular region of 
the parameter space to achieve an exact resummation of the zero modes and removing in this 
way any infrared divergence. As a result of this resummation the order parameter of the chiral 
phase transition, the chiral condensate, vanishes in the chiral limit. This signals the expected 
recovery of chiral symmetry when the quark mass m q vanishes in a finite volume, i.e. the absence 
of spontaneous chiral symmetry. When studying the cutoff effects with Wilson-type fermions in 
the framework of #PT the power counting of the £-regime in the continuum 

i = 0(£), i = 0(e), Ml = 0(e 4 ) [or m q = 0(e 4 )], (2.1) 

has to be augmented to include the lattice spacing a. This is usually done connecting the power 
counting of the quark mass m q with the lattice spacing times the appropriate powers of the QCD 

'in the text with £ we indicate, unless appearing in renormalization group invariant combinations, the chiral con- 
densate renormalized in the MS scheme and a scale of 2 GeV. 
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scale A to restore the proper dimensions. Typically in the ^-regime there are two different power 
countings depending on how the quark mass and the lattice spacing are related: the GSM regime [ |l0| ] 
and the Aoki [ JTT| ] or large cutoff effects (LCE) regime. In the so-called GSM regime ra q ~ aA 2 
which implies a = 0(e 4 ). In this regime the leading order (LO) cutoff effects can be reabsorbed 
in the definition of the quark mass, implying no cutoff effects up to NLO order corrections (this is 
indeed true also in the /^-regime). At NLO the chiral Lagrangian describing the dynamics of the 
Goldstone field U (x) is given by 



(4) 



Wx 



ifj 4) + aWTridpU* dpU)Tr(U + W 
a 2 W' [Tr(U + U^)] 2 - 2aB H'Tr{^(' + 



2aB WTr(^ n U + U^')Tr(U + t/ 1 ) + 



(2.2) 



where * s tne con ti nuum NLO chiral Lagrangian, W, W and W' are LECs parametrizing cutoff 
effects, and is the shifted mass matrix which reabsorbes the LO 0(a) cutoff effects. Given 
our particular choice of the power counting it is easy to see that all the corrections terms to the 
continuum Lagrangian are of NNLO [@, ^]. This implies that Wilson fermions are "effectively" 
free from discretization errors up to NNLO corrections. This result is indipendent whether we use 
a clover term in the fermion action or not and whether we use twisted mass fermions or not. 

The Aoki regime is defined by m q ~ a 2 A 3 which implies a = 0(e 2 ). The LO chiral Lagrangian 
in the Aoki regime is 

jgg = ^-Tr [dpU%U] - -Tr [^ n U + Jt'V*\ - a 2 W' [Tr (U + U*)] 2 . 

It contains already at LO 0(a 2 ) cutoff effects that cannot be reabsorbed in the definition of the 
quark mass. It is well known [12, |l3|] that these LO cutoff effect in infinite volume change the 
vacuum structure of the effective theory leading to two possible scenarios for the chiral phase 
diagram [IP]]. In the Aoki scenario the pattern of spontaneous symmetry breaking changes the 
continuum SU (2)z, x SU (2)r — > SU (2)y into SU (2)y — > U(l) signalling the spontaneous breaking 
of flavour (and parity) symmetry. In the Sharpe-Singleton scenario there is no phase transition in 
the chiral limit and the Goldstone bosons remain massive. It is quite clear that in both scenarios the 
physics of the zero-modes is quite different from the one of the continuum or of the GSM regime. 
More work is needed in order to understand Wilson-type fermions in the deep chiral regime. 

The £-regime gives us the possibility to have a transition region between the two regime in 
which the cutoff effects appear at NLO [§,§]. To understand this regime from a power counting 
point of view we can decide to set a = 0(e 3 ). The LO action contains, as in the continuum, the 



mass term 



S^ ' I and the kinetic term for the non-zero-modes. First corrections due to a finite 



lattice spacing appear at NLO. The partition function of the effective theory at NLO can be written 

as 



2? = yK 
where 2 
1 % [Uq] 



'[£/()] e 



- s 



J 0)1 



(Eeff) 



x Z n [U { 







,JV \ l+W'a 2 V 



TrlUo + U, 



AO) 



2 2aWLV 
+ - Tr 



-| / r/ 4 .vTr 



(2.3) 



F 2 



Tr 



(2.4) 



S e ff contains the NLO correction to E coming from the one loop non-zero modes contribution. 
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It is clear from eqs. (|2.3[) (Of) that the discretization errors of 0(a 2 ) are of NLO while the dis 



cretization errors of 0(am q ) are of higher order being of 0(£ 3 ). With this partition function it 
is straightforward to introduce appropriate sources and compute two-point functions. We refer 
to [[| ^] for details on the computation. It turns out that the relative 0(a 2 ) cutoff effects cor- 
rections, over a wide range of values for m^LV, are at most of few percent. Another result of 
our analysis []8| |9|] is that there are appropriate linear combinations of correlation functions like 
Cs(^o)/4 + 3Cp(^o) or Caa(^o) + Cw(*o) which are free from 0(a 2 ) effects and have leading cut- 
off effects of 0(am q ). We can summarize the results of our analysis generalizing the e expansion 
using Wilson fermions in the following way. In the Aoki regime there could be large cutoff effects 
and more work is needed to completely understand the interplay between quark mass and lattice 
spacing effects. In the GSM regime Wilson fermions have no cutoff effects up to NNLO. In the 
transition region between the two regimes we have a tool to analyze cutoff effects and we have 
computed correlators including 0(a 2 ) and 0(am q ). The proper power counting has been identified 
as a = 0(e 3 ) and two-point functions have been computed up to relative 0(£ 3 ) corrections. We 
have noticed that certain linear combinations have no discretization errors up to NNLO. This an- 
alytical effort obviously has to be combined with numerical simulations which can tell us which 
regime has been properly matched. We are currently extending this computation for Wilson-twisted 
mass fermions Ji5[|. 

3. Numerical results with Wilson twisted mass fermions 

To perform simulations in the £-regime with Wilson-type fermions we need specific algorithmic 
improvements. The main ingredient for these improvements is reweighting. We first proposed to 



use the PHMC algorithm [16] combined with exact reweighiting to include in an exact way the low 



modes of the Wilson twisted mass operator 060. Different but somehow related techinques have 



been later proposed based on stochastic reweighting in the standard mass [17] or in the twisted 
mass ]T8|]. The main goal of reweighting in the £ -regime is to ensure a better sampling of the con- 
figuration space and to avoid instabilities issues with HMC-like algorithms. The usage of twisted 
mass is particularly beneficial because it provides a sharp infrared cutoff for the spectrum of the 
lattice operator. We have performed simulations with a tree-level improved gauge action [ |l^ ] and 
Nf = 2 Wilson twisted mass fermions [^-22]. We summarize in tab. |]the parameters of the nu- 
merical simulations we are currently performing. To determine the LECs we compute two-point 
functions and we compare the Euclidean time dependence of the correlation functions with the 
time dependence predicted by #PT. In fig. [I] we show the numerical results, for the ensemble F 2 , 
for the charged pseudoscalar density two-point function and the fit results obtained using the NLO 



formula in the continuum [123 , 24] 



3 / \X 2 (z)\T \ ry nrv (, N lZlr(^ 

CP(X0) = ^ \xj£j + 3*« J L? hl{X0/T) f ' ^ = 2 ^ LV ^ ~ A^ G(0) J 

(see ref. []|] for unexplained notations) in the fit range 10 < xq/ci < 38. It turns out that the fit 
results are very stable if we change the number of data points included in the fit. We have repeated 
this analysis for all the simulation points of tab. [I] and the results for the chiral condensate £ and the 
pseudoscalar decay constant F have been collected in fig. ||. This figure summarizes the results of 
simulations at 2 values of the lattice spacings, 2 different physical volumes and several values of the 
quark masses, thus allowing us to attempt a first understanding of the systematic errors. Together 
with the numerical data we also plot two vertical dashed (j3 = 3.9) and dotted (/$ = 4.05) lines, 
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lattice 


P 


Lattice 


a [fm] 


ajj. q 


Ii 


3.9 


16 3 x 32 


0.079 


0.0005 


h 


3.9 


16 3 x 32 


0.079 


0.00075 


h 


3.9 


16 3 x 32 


0.079 


0.001 


F, 


4.05 


24 3 x 48 


0.063 


0.00039 


F 2 


4.05 


24 3 x 48 


0.063 


0.00078 


F 3 


4.05 


20 3 x 40 


0.063 


0.00039 



Table 1: Simulation parameters of the runs performed 
where j3 = 6/g$, the lattice spacing a J25[], /x q is the 
twisted mass parameter and iVtraj is the number of tra- 
jectories (with unity trajectory length) excluding the 
thermalization process. 



Group 


rd 1 / 3 


r F 


This work 
ETMC [f2|] 

HHS @ 
JLQCD [27] 
JLQCD [28] 


0.595(12) 
0.574(28) 
0.617(15) 
0.596(10) 
0.624(17) (27) 


0.224(12) 
0.183(8) 
0.224(10) 
0.217(14) 



Table 2: Table comparing the results presented in 
this proceedings with recent results obtained in the e- 
regime using N{ = 2 overlap fermions and clover-type 
fermions. The result obtained in the /^-regime with 
Nf = 2 Wilson twisted mass fermions by ETMC is also 
added for comparison. 



2500 
3500 
3135 
1755 
2316 
2500 



for the two different lattice spacings, indicating the value of /i q Xy where we enter the so called 
Aoki regime (see sect. |||). This value depends on the value of an unknown low energy constants, 
called C2 °= —W which parametrizes the 0(a 2 ) effects. The two lines are the bound for /XqZV for 
two different lattice spacings given the indicative, but plausible, value \ci\ = (400MeV) 4 . The 
figure can be interpreted in the following way. The two different lattice spacings and same physical 
volumes (filled simbols) agree with themselves if we exclude the most chiral point at /3 =3.9 which 
is beyond the dashed line we have drawn. This might indicate that cutoff effects might be small if 
we keep the value of jU q iy larger than the 0(a 2 W'V). The second effect is visible when we change 
the physical volume. There seems to be a discrepancy between the two volume simulated both for 
the chiral condensate and the decay constant. This discrepancy might indicate the inadequacy of a 
NLO fit when the volume is not large enough. On the other side there is a good consistency in the 
LECs if we change the value of the mass at the largest volume available for both lattice spacings 
(excluding the most chiral point already dicussed). If we take the finest lattice spacing and the 
largest volume as our best estimate of the LECs we obtain the preliminary results 

r £ 1/3 = 0.595(12), r F = 0.224(12) 

In tab. || we compare our determination for the LECs with other determinations obtained in the 
£ -regime [^-^|] and with the determination obtained by ETMC in the ^-regime. We observe a 
very good agreement with all the determinations for the chiral condensate. For the pseudoscalar 
decay constant we observe a good agreement among all the determinations in the £ -regime, while 



there is some tension with the ETMC determination in the />regime [25] 



4. Conclusions 

Probing the £ -regime with Wilson-type fermions is possible. We have now an analytical tool which 
allows us to study the combined lattice spacing, volume and mass dependence of correlation func- 
tions in the £-regime, and we have computed two-point functions at NLO. We have performed 
several simulations in the £-regime with Wilson twisted mass, using a PHMC algorithm combined 
with exact reweighting. This allows in principle the determination of the LECs £ and F without 
contaminations from chiral logs. Systematic uncertainties can now be addressed because of the 
rather large set of simulations we have performed, indicating an 8% of systematic errors coming 
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NLO chiral fit 



^ IE l = 264(10)(10) McV J 
\ 



^^^^^ 



0.72 



0.68 - r £ 
0.64 - 
0.6 - 
0.56 - , | 



I ( 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



10 20 30 40 



a ^ = 3.9 L/a= 16 
• ^ = 4.05 L/a = 20 
0.2 - ° p = 4.05 L/a = 24 



|jXV - 



Figure 1: Euclidean time dependence of the pseu- 
doscalar density two-point function and the fit results 
obtained using the NLO formula in the continuum. 
The fit curve becomes dashed outside the fit range. 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

Figure 2: Plot of ro|£| 1//3 and rQF for all the simu- 
lations point we have at the moment as a function of 



from a too small volume and a 4 — 10% of systematic errors coming from the finite cutoff. To en- 
sure a reliable understanding of the systematic uncertainties and thus remove them we will enlarge 
the set of simulation points. We will also try to enlarge the set of physical quantities to compute 
to fully take advantage of the expensive dynamical simulations perfomed so far. From the ana- 
lytical side we are curently extending the NLO computation performed with Wilson fermions to 
Wilson twisted mass [|l5|]. We see no reason why this computation could not be done with stag- 
gered fermions, and moreover it could be a good corner of the parameter space where to test the 
rooting approach. 
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